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Abstract 

We solve the mean-field (/> 4 model in an external magnetic field in the microcanonical 
ensemble using two different methods. The first one is based on Rugh's microcanon- 
ical formalism and leads to express macroscopic observables, such as temperature, 
specific heat, magnetization and susceptibility, as time averages of convenient func- 
tions of the phase-space. The approach is applicable for any finite number of particles 
N. The second method uses large deviation techniques and allows us to derive ex- 
plicit expressions for microcanonical entropy and for macroscopic observables in the 
N — » oo limit. Assuming ergodicity, we evaluate time averages in molecular dynam- 
ics simulations and, using Rugh's approach, we determine the value of macroscopic 
observables at finite N. These averages are affected by a slow time evolution, often 
observed in systems with long-range interactions. We then show how the finite N 
time averages of macroscopic observables converge to their corresponding N — > oo 
values as iV is increased. As expected, finite size effects scale as N^ 1 . 
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1 Introduction 



The microcanonical ensemble can be considered, in several respects, as the ba- 
sic ensemble for the description of the statistical behavior of physical systems. 
However, the canonical ensemble is often more amenable to exact or approx- 
imate analytical treatments and, therefore, several calculations are usually 
performed in its framework. On the other hand, there are some reasons of prin- 
ciple, and also of practical importance, for which the microcanonical ensemble 
must be given a prominent role [1] . Indeed, there is an increasing evidence that 
ensembles can be nonequivalent for systems with long-range interactions [2] 
or for finite (small) systems, when surface effects are relevant [1,3]. In the 
mathematical physics literature, this question has been first discussed in the 
pioneering paper by Hertel and Thirring [4] (a recent publication where several 
references to these results can be found is Ref. [5]). An alternative approach, 
in which a parallel study of the distributions of microstates and macrostates 
is performed, has been proposed in Ref. [6]. This approach is based on large 
deviation techniques [7]. 

On the other hand, it has been shown that inequivalence between microcanon- 
ical and canonical ensembles can be associated to the presence of first order 
phase transitions in the canonical ensemble [8], and that there are regions in 
the phase diagram of the system that can be reached only within the frame- 
work of the microcanonical ensemble [9,10,11]. This means that in such cases, 
preparing the system at either fixed energy or temperature determines dif- 
ferent values of the macroscopic observables. It has been correctly pointed 
out that this implies the impossibility of defining a unique equation of state 
for the system [3]. Moreover, in systems with long-range forces or in small 
systems, when the interaction between macroscopic subsystems cannot be re- 
duced to a surface effect, the very concept of a heat bath can be seriuously 
questioned [12]. 

The physical consequences of ensemble inequivalence are of paramount impor- 
tance. Perhaps, the most striking example is the existence of negative specific 
heat in the microcanonical ensemble for self-gravitating systems [13] . This re- 
flects the real physical situation, and, in the canonical ensemble, negative spe- 
cific heat is simply impossible. Experimental observations of negative specific 
heat for atomic and molecular clusters have also been recently reported [14,15]. 

The above arguments, which are all consequences of ensemble inequivalence, 
justify the necessity to perform, at least in some cases, microcanonical ensem- 
ble calculations. The operative problem then arises on how this can be actually 
done. This question can be faced both in the thermodynamic limit and at a 
finite number N of degrees of freedom. The study of finite iV is relevant in 
two respects: i) in Molecular Dynamics (MD) simulations one would like to 
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estimate the approach to the N — > oo limit in numerical experiments, ii) finite 
(small) systems, as just noted, can be interesting by themselves. 

In this paper we discuss two different approaches to the study of a physical 
system within the framework of the microcanonical ensemble. We will consider 
the mean-field 4 model in an external magnetic field. We will here restrict 
to the description of the solution methods and to the implementation of the 
expressions we obtain for macroscopic observables in MD simulations. We will 
not discuss issues directly related to ensemble inequivalence, which will be the 
subject of a forthcoming publication [16]. 

In Section 2 we review two different approaches to the calculation of micro- 
canonical averages. The first one is based on the microcanonical formalism 
introduced by Rugh [17,18]. The second [19] relies on large deviations tech- 
niques, which are particularly suited for mean-field like systems. In Section 3 
we present and briefly discuss the 4 model. In Section 4 we show how the two 
approaches can be applied to the model. In Subsection 4.1 we derive the finite 
N formulas for some observables in Rugh's microcanonical formalism, while in 
Subsection 4.2 large deviations techniques are used to compute the entropy of 
the system, and, from it, the expressions of the thermodynamic observables. 
We then show that the large N limit of the microcanonical formalism gives 
the same results as the large deviation approach. In Section 5 we present the 
results of MD simulations of our model, with the aim of illustrating the ap- 
proach to the N — > oo limit. Section 6 is devoted to some conclusions and 
perspectives. 



2 Two ways of obtaining microcanonical averages 

In the computation of the entropy and of the averages of observables, the 
problem at hand is the calculation of integrals of functions of the phase space, 
restricted to a given energy shell, and of the dependence of these integrals on 
the energy. Here we consider two procedures that are somewhat complemen- 
tary. They have been recently introduced in the literature (Refs. [17,18] for 
the finite N approach and Refs. [19,20] for the large deviation one), and our 
aim is here to make a comparison between them, both as far as general aspects 
are concerned, and for the actual implementation in a specific case. 

Rugh's microcanonical formalism [17,18] is in principle more general, since 
it is formulated for general Hamiltonians and for any number of degrees of 
freedom N. It does not allow to obtain the entropy itself, but it gives access 
to averages of the observables and to derivatives of the entropy and of the 
averages with respect to the energy and to the parameters of the Hamiltonian. 
These derivatives are in turn expressed as averages of suitable observables, 
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which must be determined otherwise, e.g. numerically. On the other hand, the 
large deviation method [19,20] is mostly applicable to systems with long-range 
interactions. Expressions are obtained in the thermodynamic limit, although 
finite N corrections can in principle be derived as series in 1/N. Its advantage 
is the possibility to actually compute the entropy itself, besides its derivatives, 
and to obtain analytically features of the N — > oo limit, like phase transitions. 



2.1 Rugh's microcanonical averages 



In this Subsection we review the basic expressions for the averages of observ- 
ables in the microcanonical ensemble, as given in Refs. [17,18]. We restrict 
to the case in which the total energy is the only integral of motion, and we 
consider a Hamiltonian that can depend on a parameter A. 

In units where the Boltzmann constant is equal to 1, the microcanonical vol- 
ume is given by: 

Q(E, A) = e 5(£;A) = jdT5(E- H(T; A)) , (1) 



where T denotes the phase space of the system, and A the parameter the 
Hamiltonian depends on. The average of an observable A(T) is: 

/ a(\\\ SdT6(E-H(T;X))A(T;X) 



It is possible to express in turn the derivatives with respect to E and A of 
{A(X)) E ,x and of S(E;X) (and thus also the temperature) as averages. It is 
useful to introduce the notation 

/i(A(A); E,X) = JdT5(E- H(T; A)) A(T; A). (3) 



The derivatives of the function fi with respect to E and A can be expressed 
using any vector Y in T space such that Z^l^f^- = Y • 'VH = 1 (where Xi 
are the phase space variables and the gradient operator V acts on T). Using 
the property 



E Yi£-6(E - H(T; A)) = -8\E - H(T; A)) £ 

= -5\E-H(T;X)), (4) 



one gets 
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— pl(A(X); E,\) = J dT5'(E - H(T; X))A(T, A) 



dr 



£^(£-#(r;A)) 



= -/ 

= J dT5(E - H(T; A)) 



A(F;X) 



d 



= H(V-(AY);E,X), 
where we have integrated by parts. Moreover, 



(5) 



—v(A(\y,E,\) = Jar 



-5'(E-H(T;X)) 



dH(T; A) 



+S(E-H(T;X)) 



dA(T; A) 



dX 



= J dT5(E - H(T; A)) 



dA(T; A) 



-f (v ■ ( a ~^AY)\ E,\)+v(^, E, x) . (6) 



The entropy S(E; A) and the average {A(X)) E a are expressed, through the use 
of (3), by ' 



S(E; A) = ln/i(l;£,A) 
fi(A(X);E,X) 
(A(A)kA= n(l;E,\) ■ 

From Eqs. (5), (6) and (7) it follows that 



|WA)^-(V.(fY)), A 
^(A(X)) E , X = (V • (AY)) E>X - ¥ ^ ) (A)e,x 



(7) 



(8) 
(9) 
(10) 



d_ 

OX 



+ 



"(V 

,dA(X). 
v <9A ' 



OH 
~dX 



AY)) E , x +(V-(—Y)) E , x (A) E 



IE,\ 



(11) 
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In the simulations, ensemble averages are computed through time averages. 
This is allowed when the ergodic hypothesis is satisfied. However, it can be 
argued that, from the practical point of view, this is also allowed in the pres- 
ence of metastable states. In fact, in such cases the phase-space trajectory is 
confined for a long time (longer than the duration of a simulation) in a subset 
of the constant energy surface. Time averages computed along this trajectory 
can be interpreted as representative of ensemble averages in which another in- 
tegral of motion, beyond the energy, is present, i.e. the characteristic function 
of the subset where the motion is limited. 



2. 2 Entropy from large deviations 



Being our system of the mean-field type (see next Section), its micro canonical 
entropy, given by Eq. (1), can be also obtained using large deviation tech- 
niques [19]. These techniques are more generally suitable for systems with 
long-range interactions [2]. It is interesting to derive ensemble averages also 
in this framework. 

Let us just begin with a brief illustration of how entropy is computed using 
large deviation techniques. We suppose that the Hamiltonian of the system 
with N degrees of freedom can be written in the following mean-field like form: 

n 

H = Nj29k(M k ), (12) 
k=i 



with 



1 



A? 



M k = -J2h(<HiPi)i (13) 

i=l 



and where f k ,g k in (12) and (13) are smooth functions. Leaving implicit the 
dependence on quantities other than the energy, the microcanonical volume, 
Eq. (1), can be transformed as: 



e s(Ne ) = J dr5 (H(r)-E) = Jdr8(Np^g k (M k (r))-Ne^ 
= jdT [U n k=1 dm k 5(NM k (T) - Nm k )} 5 
Using the inverse Laplace transform of the 5 function we get: 



N \J29k(m k ) - e 
\k=i / 



(14) 
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e S(Ne) = J dT [U n =idmkdzk] ( JL\ l S 



\27ri 



N Eff*W ~~ e 
. \fc=i j 



exp 



2 (Nz k M k (F) - Nz k m k ) 

.k=l 

N 



= J [U n k=1 dm k dz k ] (— ) 5 



N Yl 9k(m k ) - e 



\k=i 



exp <^ -JV 



J2z k m k ) - ln(e^^i 2fe/fc )i 



L \fc=l 



where (-)i is the one particle phase-space (q,p) average given by: 

(•>! = / dqdp(-) , 



(15) 



(16) 



which is of course defined only for integrable functions. In Eq. (15) the integra- 
tion path for the variables z k is a line parallel to the imaginary axis and with 
a positive real part. The integrals in z k are computed, for iV — > oo, using the 
saddle point method. The relevant stationary point of the function in square 
brackets in the exponent in (15) must lie on the real z k axes. Otherwise one 
would find an unphysical oscillatory behaviour of the integral. The minimum 
along a line parallel to the imaginary axis is a maximum along the real axis; 
therefore we get: 



N \J29k(m k ) - e 



\k=i 



exp [-NI (mi, . . . ,m n )] , 



where: 



I (mi, . . . ,m n ) = sup 

Zl,...,Z k 



1 a \ 

J2z k m k ) -ln(e^i 2fc/fc )i 



\k=i 



(17) 



(18) 



with the sup taken for real values of the z k s. Finally, in the iV — > oo limit we 
have: 



s(e) = lim ^S(Ne) = -1(e) , 



(19) 



where: 



/(e) = inf I(m 1 ,...,m n ) = 

mi,...,m n ,229k(mk)=e 
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sup 

<,^9k{m k )=e 



[-/(mi,...,m„)] . 



(20) 



mi,...,m„, 



If the last inf is taken only on some of the m&, then we obtain the entropy as 
a function of e and of the remaining m^'s, which are then considered as given 
quantities. 

Later, in Section 4, we will implement both procedures in our system. Since 
expressions computed using large deviation techniques describe the system in 
the N — > oo limit, the comparison will be performed with formulas obtained 
with Rugh's microcanonical formalism taken in the same limit. 

The important point to be noted is the following. Rugh's formalism is more 
general, since it is valid, in principle, for any system and for any size N . It 
provides expressions for the average of observables and for the derivatives of 
these averages, that afterwards can be effectively computed in simulations. 
However, the microcanonical volume, and the entropy itself, can not be evalu- 
ated. On the other hand, large deviation techniques are more easily applied to 
systems with long-range interactions, but one can obtain an explicit expres- 
sion of the entropy, which gives access to the full knowledge of thermodynamic 
properties, including phase transitions. 



3 The 4 model in an external field 

The Hamiltonian of our system is given by: 



where the qi variables can be thought as describing the position of the i-th 
particle on a line, and the p^s are the conjugate momenta (the mass is unitary). 
Each particle is subject to both a local potential and to an infinite range one, 
expressed by the all-to-all coupling; h is the external magnetic field and 9 is 
a free parameter. Indeed, it can be shown that all parameters appearing in a 
potential energy of this form can be absorbed in 9 by a convenient change of 
variables. In the MD simulations we have chosen the value 6 — 1/2, which is 
in the range that ensures the presence of an effective double well potential in 
the low energy phase. The Langevin dynamics with the force coming from the 
potential energy of model (21) has been solved in Ref. [21]. A careful study of 
the different dynamical regimes for the relaxation to equilibrium and a first 
characterization of the second order phase transition has been performed in 
the same paper. A further study of the canonical ensemble solution of this 
model for h = has been performed in Ref. [22], where it has been shown 
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that the system displays, for 9 > 0, a second order ferromagnetic transition. 
The magnetization defined by 

1 N 

M 1 = -]>> (22) 
JV i=i 



vanishes continuously at a critical temperature, which can be determined nu- 
merically by solving an implicit equation. In Ref. [22] the authors have shown 
that for 9 = 1/2 the critical temperature is T c ~ 0.264. For this particular 
value of 9 it can be shown that, curiously, the corresponding critical energy 
density e c is equal to T c /2, thus e c ~ 0.132. 



Hamiltonian (21) can be cast in the form (12): 



H = N 



P - (hM 1 + °-Ml\ - \{l - 9)M 2 + ]m a 



(23) 



using the n = 4 mean- fields (see Eq. (13)) 

1 N 



i=i 



N 



iV i=l 



(24) 



with s = 1,2,4. 



4 Implementation of the general expressions 



In this Section we will implement the general expressions introduced in Sec- 
tion 2, discussing the solution of the mean-field 4 model in the microcanonical 
ensemble. 



4-1 Finite N microcanonical averages 



Following Ref. [17], we take for Y the vector given by: 

Y= jy 2 (pi,...,Piv,0, ...,0) = ^(pi,---,Pjv,0, ...,0), (25) 



with non vanishing components only in correspondence of the momenta. The 
external magnetic field h plays the role of the parameter A. We are interested 
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in the expressions for the temperature, the specific heat, the average magne- 
tization and the magnetic susceptibility. Computing V • Y, the temperature 
follows from Eq. (8): 

1 .N-2, . . 

-)E,h- (26) 



T(E,h) v 2K 



From this, and using (9), the inverse specific heat is given by: 



1 



C(E, h) 



m T ^ = 

l-T 2 (E,h)( 



d 



1 



~ T2{E ' h) dET(E,h) 
(N-2)(N -A) 



(27) 



AK 2 



>E,h 



= 1 



(N 



4) (l/K*) Eth 



(TV -2) (l/K)l :h 



We point out that C is the extensive specific heat and thus Eq. (27) is expected 
to scale as N" 1 . The average magnetization {M-\)e,Hi with Mi given by (22), 
is denoted by m(E, h). Then the magnetic susceptibility at constant energy is 
computed through (10), noting that dH/dh = —NMi, and it is expressed by: 



X (E,h) = ^m(E,h) = N(N-2) 



Mi 

(^)E, h m(E,h) 



■ (28) 



4-2 Microcanonical solution of the (ft 4 model 



In this Subsection we will present the solution of the mean-field 4 model using 
large deviations techniques. In particular, we will derive explicit expressions 
for the microcanonical entropy and its derivatives in the iV — > oo limit. We 
will apply the formalism introduced in Subsection 2.2. 

The one particle average in the rightmost side of Eq. (15) takes the form 



ip(z p , zi, z 2 , z 4 ) = J dqdp exp 
\J -z p I ^ q 



z p p 2 + Zl q + z 2 q 2 + z,q 4 



exp 



z x q + z 2 q 2 + z 4 q 4 



(29) 



where convergence requires that z p and z^ be negative. The function i/)(z p , Zi,z 2 , Z4) 
is the so-called generating function in large deviation theory. 



Therefore, the extremal problem defined by (18) is: 



I(u,Q 1 ,Q 2 ,Q A ) = sup [zpU + Z]Q\ + z 2 Q 2 + z a Qa - hiip(z p , z x , z 2 , z A )\ , 

Zp,Zl ,22,Z4 
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(30) 



where the symbols u and Qi, Q 2 and stand for the in (15). The vari- 
able z p separates from the others and, after straightforward calculation, one 
obtains: 



I(u, Qi, Q 2 , Qi) = ~ ln(27re) - ^ lnu + z 1 (Q 1 , Q 2 , Q^Qi + 

+ Z 2 {Q 1 ,Q 2 ,Q i )Q 2 ^Z A (Qi,Q 2 ,Qi)Q A -h\G{z 1 ,Z 2 ,Z A ) , (31) 

where the function G is given by: 

G(zi, z 2 , z A ) — dq exp z x q + z 2 q 2 + z A q 4 , (32) 



and the functions z s (Qi, Q 2 , Q4) are the solution of: 

Qs = ~q~q^~^'( Zi > z<2 > z ^ = ~&( Zu Z<2 ' Z ^ 



(33) 



for s = 1, 2, 4. Here G s is expressed by: 



G s (z 1 ,z 2 ,z 4 ) = J dqq s exp 



z x q + z 2 q 2 + z^q 4 



(34) 



In the actual computation one has to check that the stationary points found 
from (33) correspond to the proper extrema. Finally, using (19) and (20), the 
entropy is: 



s(e,h)— sup {-In 
Qi,Q 2 ,Q4 



2e + (2hQ 1 + 6Q\) + (1 - 9)Q 2 - -Q 4 



+ 



+ In G(z ± , z 2 , z A ) - Z1Q1 - z 2 Q 2 - z±Q^ . 



(35) 



Here we have used the equation of the total energy to eliminate the variable 
u. Using Eq. (33), that defines the functions ~z s {Q\iQiiQi)i the stationary 
points that solve the variational problem (35) are: 



0Qi + h 

= Zi 



u 

1-9 



= z 2 



~Tu~ z ' 



(36) 
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where, for brevity, the expression in square brackets in (35) has again been 
denoted as u. The solution for Q s (and u) can be denoted with Q s (and u). 
It can be seen that, as expected, z 4 is a negative quantity, since the kinetic 
energy per particle |TZ is positive. Also, one has that z 2 = — 2(1 — 0)z 4 . 

The first derivatives of the entropy can now be explicitly computed, using well 
known thermodynamic relations. Taking into account the stationary point in 
Eqs. (33) and (36), one gets: 



ds(e, h) 


_ 1 


1 


de 


= T 


u 


ds(e, h) 


m 


= Qi 


dh 


= T 


u 



(37) 



The first equation is exactly the same as Eq. (26), which is obtained using 
Rugh's formalism, once it is taken in the thermodynamic limit. The second 
of the Eqs. (37) is practically a tautology, since it can be shown that the Q s 
are simply the averages of the quantities M s defined in Eq. (24). We also note 
that the functions ~z 2 and z 4 can be expressed as ~z 2 = (1 — 0)^ and Z4 = — ^ • 

Before considering the second derivatives, we show how to obtain the critical 
point. We consider the system with h — 0. Then, Eqs. (33) and (36) have 
always the solution z± — Q 1 — 0. However, for an energy density e smaller 
than a critical energy e c , the relevant solution is one with Q l > 0, although 
we do not prove it here. We only show how to obtain the value of e c and the 
corresponding critical temperature T c . 

From Eq. (36) one gets the relations z 4 = — z 2 /(2 — 29) (already mentioned) 
and z± = 1QQ{z2j (1 — 0) (for h — 0). One can insert these relations in Eq. (33). 
Near the critical point, Q l is small; we can therefore write Eq. (33) for s — 1, 
performing a power series expansion in Q 1 up to first order. Calling z 2c the 
critical value of z 2 , we easily obtain the following equation for z 2c : 



j^ 2c / dg g 2 exp |z 2c 


V 2 2 { Le)A} 


fdq exp {z 2c 


a 2 1 o 4 " 

y 2(1-0) y _ 


} 



= 1 , (38) 



which can be easily solved numerically. This same equation also tells that: 

= (39) 

It can be also easily shown that: 

Q 4c = Q 2c . (40) 
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From the second of Eqs. (36) and from the Hamiltonian, the critical temper- 
ature and the critical energy are found to be : 



T r = u r 



1 - 9 

2Z2c 

;i-fl)(2fl-j) 



(41) 



Let us now consider the second derivatives of the entropy and, to be definite, 
let us consider the specific heat computed from the entropy (35). We do 
not show here the cumbersome expression obtained differentiating twice the 
entropy (35) with respect to e. However, it can be easily guessed that several 
terms containing the coordinates, and not only the momenta, will appear. 
This is in contrast with Eq. (27) derived from Rugh's formalism, where only 
the kinetic energy appears. Of course the results must be equivalent, but 
the connection between kinetic and potential energy in the microcanonical 
ensemble (absent in the canonical ensemble) hinders the equivalence. However, 
an equation analogous to (27), when this is taken in the thermodynamic limit, 
can be still obtained if one performs the derivatives before performing the 
integrations leading to Eq. (35). Let us focus on the 5 function in the rightmost 
hand side of (15). For our system its argument is: 



N 



e 



-u 



Ql 



0)Q2 + \Qa 



hQ l 



(42) 



Thus, the derivative with respect to e of the 5 function is equal to —2 times 
the derivative with respect to u. Integrating by parts in u we then find that: 



d 2 S{Ne) 
d{Ne) 2 



(4) - ( z P y 



(43) 



with the microcanonical average (•) computed through the microcanonical 
formula (15). The saddle point evaluation of (43) has to be performed up 
to order 0(N~ 1 ), since the terms of 0(1) cancel out. The calculation, after 
several steps, leads to: 



d 2 S(Ne) 
d(Ne) 2 



+ 



((5u) 2 ) 



N{u) 2 (u) 



(44) 



where u, as before, is twice the kinetic energy, and ((Su) 2 ) denotes the variance 
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of u. Therefore the inverse specific heat is: 




g d 2 S(Ne) = 2_ _ ((^) 2 ) 
d(Ne) 2 N (u) 2 



(45) 



In (44) and (45) the two terms are both of order iV -1 , in spite of their ap- 
pearence. It can be easily shown [17] that the N — > oo limit of (27) leads to 
exactly the same expression. 

The same kind of procedure can be followed for the magnetic susceptibility. 
Now we can use the fact that the derivative with respect to h of the 5 function 
in the rightmost hand side of (15) is equal to — 2Q 1 times the derivative with 
respect to u. Integrating again by parts we get: 



which is the same as the thermodynamic limit of (28), reminding that Qi 
corresponds to the dynamical variable Mi given in Eq. (22). 



5 Molecular dynamics simulations 

In this Section we show some results of MD simulations of our system. In 
particular, we present numerical data on the microcanonical temperature, the 
specific heat, the magnetization and the magnetic susceptibility. 

The simulations have been performed with N = 25, 50, 100, 200 particles, and 
the value 9 = 1/2 has been selected. We have considered here only energies 
above the critical energy e c ~ 0.132 (see Section 3), since in the simulations 
presented in this work we only want to show that the microcanonical formalism 
can be easily applied, in practical numerical computations. Below the critical 
energy there can be metastable states [22], that we do not consider here. In 
a work in preparation, devoted to nonequivalence of statistical ensembles, we 
will also present the results of numerical simulations performed below the 
critical energy [16]. 

We have found that the convergence of the averages can be significantly im- 
proved if, during the simulation, the velocities of the particles are reshuffled 



X = -2N (Qjz p ) - m(Q lZp ) 



(46) 



The leading order in N of this expression gives: 




(47) 
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at given time intervals. This is particularly true for the quantities related to 
the second derivatives of the entropy, i.e. the specific heat and the magnetic 
susceptibility. The faster convergence is presumably due to a more efficient 
spanning of the phase space of the system. 

In Fig. 1 we show the results for the energy density e = 0.25. The simulation 
time is 2 ■ 10 5 . The figure contains four plots in log- log scale, in which we 
show, as a function of the number of particles, the difference between the av- 
erage value of four observables obtained in the finite N simulations and their 
corresponding N — > oo limit computed through the large deviation method 
(indicated by the subscript oo). The four observables are: the microcanonical 
temperature T (26), the magnetization m = {Mi)e,h (22), the inverse specific 
heat 1/c (27) and the magnetic susceptibility x (28). Actually, for the mag- 
netization, the plotted quantity is the N — > oo limit value minus the finite 
iV value, which is smaller. Furthermore, for the inverse specific heat we have 
considered the intensive quantities that are obtained multiplying (27) and (45) 
by iV (and then we have used the lower case for the corresponding symbol). 

Assuming a power law decay with iV of the simulated averages towards the 
iV — > oo limit value, the exponent of this law can be estimated from a fit with 
a straight line. In the plots the lines are interpolations between points, as a 
guide to the eye. We do not show the fitting lines that can be easily obtained. 
We limit to point out that the slope of these lines is in all cases close to 
— 1. This behavior should be expected on general grounds; in fact, the finite 
iV corrections could be obtained by the successive terms in the saddle point 
evaluation related to the large deviation computation, and the first correction 
is expected in general to scale like AT -1 . 



6 Conclusions and perspectives 

In this paper we have considered in parallel two methods for the evaluation 
of microcanonical averages of macroscopic observables and we have applied 
them to the mean-field 4 model in an external magnetic field. 

We have implemented in the model the general expressions that, through the 
microcanonical formalism introduced by Rugh [17,18], allow us to compute 
an entire set of macroscopic quantities (temperature, magnetization, specific 
heat, magnetic susceptibility, etc.), performing time averages of suitable me- 
chanical observables at finite N. In particular, these macroscopic variables are 
obtained taking derivatives of the entropy with respect to the energy density 
and to other parameters of the Hamiltonian, like the external field (see e.g. 
formulas (37)). 
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Fig. 1. Microcanonical time average of temperature T (a), magnetization m(b), 
inverse specific heat 1/c (c), magnetic susceptibility x(d). In all figures we show the 
difference of the time average at finite N with respect to the N — > oo limit value 
(which is analytically computed using large deviation techniques and is indicated 
by the oo subscript) vs. the number N of particles. The energy density is e = 0.25 

Then, we have considered large deviation techniques, whose applicability relies 
on the determination of a set of macroscopic variables on which the energy 
depends in a simple way. This is generally possible for systems with long- 
range interactions, although the method is viable also for other cases [19]. 
Our model, being of the mean-field type, can be solved exactly in the micro- 
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canonical ensemble, using large deviation techniques. This leads to the analytic 
determination of the same macroscopic quantities in the thermodynamic limit. 

The purpose was here to derive the actual expressions of macroscopic observ- 
ables that stem from the application of the two different methods mentioned 
above to a model that has already been a subject of study for its interesting dy- 
namical properties [21,22]. Although these two approaches lead to expressions 
which show no evident similarity, especially for specific heat and magnetic sus- 
ceptibility, the time averages computed with Rugh's microcanonical formalism 
converge to the thermodynamic limit value computed with large deviations. 
Corrections are of order A -1 . 

We find that the specific heat and the magnetic susceptibility converge ex- 
tremely slowly in time, unless one devises methods to "help" the system ex- 
plore more efficiently phase space. The method we use consists in periodically 
reshuffling particle velocities. These observables, which correspond to higher 
order derivatives of the entropy, have extremely complex expressions, and this 
could be at the origin of the observation of such slow convergence. However, 
we believe that this slow relaxation is rather a distinctive signature of the 
long-range (mean-field) nature of the model. Indeed, slow time relaxations 
are quite common in such systems and have been already observed in other 
models [2,22,24,25], being associated to the presence of the so-called quasi- 
stationary states. 

In a future paper [16] we will concentrate on the properties of this model in 
the low energy/temperature phase. We will show the presence of ensemble 
inequivalence features and of a negative susceptibility in the microcanonical 
ensemble. 

While finishing this paper we became aware of a very recent study of the 
4 model without external field which uses large deviation techniques [23]. 
However, this work focuses on different aspects, being devoted to a detailed 
study of the non-analiticity properties of the entropy function. 
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